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Abstract —In this paper we examine convergence properties of 
massive MIMO systems with the aim of determining the number 
of antennas required for massive MIMO gains. We consider three 
characteristics of a channel matrix and study their asymptotic 
behaviour. Furthermore, we derive ZF SNR and MF SINR for a 
scenario of unequal receive powers. In our results we include 
the effects of spatial correlation. We show that the rate of 
convergence of channel metrics is much slower than that of the 
ZF/MF precoder properties. 

I. Introduction 

Multiple-Input Multiple-Output (MIMO) has been incorpo¬ 
rated into emerging wireless broadband standards such as LTE. 
In general, in time-division duplexing (TDD) mode, as we add 
more and more antennas to the base station, the performance 
in terms of data rate, enhanced reliability, improved energy 
efficiency and interference increases ID. Consequently, mas¬ 
sive MIMO is an emerging technology, where the number of 
antennas is scaled up by 1-2 orders of magnitude El-ill with 
the aim of even greater performance benefits ID. 

Motivating the surge in research activities into massive 
MIMO are the additional gains resulting from a large channel 
matrix, due to the asymptotics of random matrix theory. For 
example, the effect of small-scale fading can be averaged out 
0, beacuse random elements tend to become deterministic 
and matrix operations can be computed much easier due to 
matrices becoming well conditioned 0. Consequently as we 
let the number of antennas tend to infinity, system analysis 
becomes simplified. 

However, as a very large number of antennas may be 
impractical due to physical constraints, we would like to 
determine the point where MIMO systems begin to exhibit 
these additonal benefits. This is the focus of the paper. The 
following are the contributions: 

• We study the convergence properties of the channel 
matrix by considering three metrics defined in Sec. III. 

• We investigate precoder properties by deriving new an¬ 
alytical expressions for SNR and SINR of ZF and MF 
precoders for unequal power gains. 

• We demonstrate the difference in convergence rate of 
channel properties and that of precoder characteristics Q. 

II. System Model 

A. System Description 

We consider a massive MIMO system with M co-located 
antennas at the base station serving K single-antenna users. 

'Note: Imperfect CSI is considered in future work 


On the down-link (where TDD is assumed) the K terminals 
collectively receive a K x 1 vector 

x /= v / P/ gTs /+ w /> (!) 

where pf is the transmit SNR, [-] T represents matrix transpose, 
s/isaAfxl precoded vector of data symbols, q/, and w/ is a 
Kx 1 noise vector with independent and identically distributed 
(iid) CJ\f(0, 1) entries. The transmit power is normalized such 
that E{||sj||} = 1. The M x K channel matrix, G, is 

G = HD|, (2) 

where H is a M x K matrix which accounts for small-scale 
Rayleigh fading and transmit spatial correlation, and is 
a diagonal matrix modelling large-scale effects. The diagonal 
elements of D« are represented by a K x 1 vector, /3, with 
entries 8j representing the link gains. 

B. Power Model 

We consider cases with equal and unequal link gains. The 
equal power case models a single-user MIMO system where 
one user has K co-located antennas. This is used as a reference 
case. While not considered by other authors, the unequal 
power case models K distributed UEs where each user has a 
different link gain due to path-loss, shadowing, etc. Since we 
are considering convergence issues as the system dimension 
grows large, it is not convenient to generate the link gains 
using classic log-normal shadowing, path loss models. This 
is because the variation in the link gains will confound the 
limiting effects. Also, with random link gains, as the system 
size increases we may obtain an artificially large number of 
high gains or small gains. We propose a power model to 
counter these two problems. 

We select the 8 j values from the limiting function f3(x) = 
Ap x , where p is arbitrarily set such that 0 < p < 1, A = 
/3 max , an d 0 < x < x 0 , such that xq = 1 ° s ^ 1 "g n (^ max ' ) . This 
allows us to control the range of the link gains in the interval 
[8mm, Anax] = [Ap x °, A] and also to control the rate of decay 
of the link gains by p. Given the parameters 8 mia, /3 max , V and 
the number of users, K, the /3j values are given by 8 j = 
— 1)) which gives the K values of 81 , 82 , ■ ■ ■ ,8k 
as the values of 8(x) using K values of x evenly spread over 
[0, xo]. 

C. Convergence Model 

Many features of massive MIMO systems are driven by the 
convergence of various functions as the system size increases. 



A key issue is the convergence of jjG r G“ to I. This is exam¬ 
ined in Sec. IIII- Al and IIII-Bl Also of interest is the convergence 
of ZF and MF precoder properties to their asymptotic limits. 
Thus, in Section IIII-CI we examine the convergence of the 
SNR and SINR of ZF and MF precoders, respectively. 

In doing so we will consider the following two types of 
convergence scenarios: 

1) K is fixed, M —> oo, 

2) -jg = a~ r ,K —> oo, a is fixed. 


D. Correlation Model 

In massive MIMO, if we deploy more and more antennas 
in a fixed volume, correlation amongst antenna elements will 
increase. Hence, as we allow the array sizes to grow, we 
also allow correlations to increase. We consider the Kronecker 
correlation model, where 

H = R t 1/2 H iid , (3) 

where Hud is the M x K iid channel matrix, H is the M x I\ 
correlated channel matrix and R t is the M x A/ transmit spatial 
correlation matrix. The elements of R t , r,; ; -, are based on the 
exponential correlation matrix model Col 

Uj = P dij , \p\ < 1, (4) 

where dij is the distance between antennas i and j, p is the 
correlation decay constant and is the correlation coefficient 
between antennas i and j. We will assume the M colocated 
antennas are in a Uniform Linear Array (ULA) at the base 
station. 

III. Methodology 


In this section we focus on the second convergence scenario, 
where both K and M grow large with a = M/K remaining 
fixed 0 


A. Convergence of G T G* /M 


First, consider the entries of H 2 H*. Let H = 
[hi I 12 ... h k\, where R is an M x 1 column vector with 
entries, hij, then it is well known that 


H T H* \ hfh* E^lilM 2 X 

M ).. M M M’ 

/ ll 


(5) 


where X is X% M with 2 M complex degrees of freedom. 
Similarly 


H t H*\ hfh* h ir h* r Y 

M ),, M M ~ M' 

' L J 


( 6 ) 


Using the known convergence of -jjH v H* to I we have | [TT1 


G t G* 

M 


I) H'lT T) : 

M 


—>• Dc 


(7) 


2 Due to limitations of space all the analysis shown in Sec. lIIIl is for the iid 
case, but can be easily modified to include the effects of correlation. 


where DJg? is the limiting value of assuming it exists. The 
diagonal elements of -^G T G* are given by 

(G T G*)ji p t X 

M ~ ~M (8) 

while the off-diagonal components are given by 


(G r G%- 

M 




(9) 


Expanding X and Y in terms of the complex Gaussian 
variables, hij, using ([3 and ([6]), it is straightforward to show 
that: E(X) = M, E(F) = 0. Furthermore, Var{jj} = 
Var{jj} = jj, and thus when M — > oo the variance 
becomes zero. Hence, we have E {-j(G T G*)ij} = Pi, 
E{^(G T G*) y -} = 0, Varjc^Gl^} = pJM and 

Var |( G J 3 )ij | = yj'Pipj/M. Essentially, the speed of 
convergence of -A G 7 G * is controlled by the convergence of 
w = ^h t h* to I with the Pj values scaling the variances 
of the elements of W. 


B. Convergence Metrics 


We can further evaluate the convergence of W by examin¬ 
ing a number of well known properties of W and a deviation 
matrix E = W — I. Letting Ai, A 2 ,..., A k denote the 
eigenvalues of W, we consider the following metrics: Mean 
Absolute Deviation (MAD), A ratio and Diagonal Dominance, 
defined as 


MAD(E) = -T J2 

i=l,j=l 


E 




A ratio = 


Diagonal Dominance = 


Amax(W) 


An 


»(W)’ 


Z—/ 7— 


■3=1,3 & l Wi jl 


GO) 

( 11 ) 

( 12 ) 


These metrics will be evaluated via simulation for a number 
of system scenarios in Section ITVl 


C. SNR/SINR convergence for ZF and MF precoders 


In this section, we derive limiting expressions for the ZF 
SNR and MF SINR for the unequal power scenarios. We 
include a summary of the equal power results (9} as these are 
needed for the derivations. As in 0, we consider convergence 
scenario 2. 


1) Zero Forcing Precoding , Equal Link Gains: For a zero 
forcing precoder, the transmitted symbol vector, Sf is given 
by 

s/ = —G , (G T G’) _1 q/, (13) 

VT 

where we normalize the average power in Sf to pf via 
tr((G T G*)- 1 ) 


7 = 


K 


(14) 


























Note that the entries of G are iid in this case. The resulting 
instantaneous received SNR is given by 

Pf 


SNR = 


tr((G T G*) 


-if' 


The limit of (IT5] > as M, K —> oo, with fixed ^ = a is ||9l 


SNR —»• p/(a — 1). 


7 = 


tr((D|H T H*D|)- 1 ) 
K 


Hence, as in (fl5] >. the instantaneous SNR is given 

Pf 


SNR = 


tr((D|H T H*D|)- 1 ) 


E 


{(d|h t h*d /3 )- 1 } 


M-K' 


J_ (1 

A?=1 Bi l B 


Bi 


B J 


M-K a- 1 ’ 


SNR 


G) c 


3) Matched Filter Precoding, Equal Powers: The transmit¬ 
ted signal for a MF precoder is given by 


s, - iov 


with 


tr(G T G*^ 
~K 


7 = 


The received signal is thus given by 

x/ = v^7G T G*-^+w /l 
giving the instantaneous SINR of the <th user as 


SINR; = 


^\sfs*\ 2 


1 + jfr gfgfcgfc g; 




M 


= a fixed 


SINR, 


Pf a 

Pf + 1 


4) Matched Filter Precoding, Unequal Powers: Here, the 
power normalization factor is 


(15) 


7 = 


tr(D|H T H*D|) 
K ' 


(27) 


(16) 


2 ) Zero Forcing Precoding, Unequal Link Gains: General¬ 
izing the analysis in Section IIII-Cll to unequal link gains, we 
have the power normalization 


From (125b . the instantaneous SINR, for unequal powers, of 
the ith user can be shown to be 

^A 2 |hfh*| 2 


SINR, = 


l + ^£f=Wi&l h fh*| 2 - 


(17) 


(18) 


In order to examine the asymptotic behaviour of 
rewrite the right-hand side to give 

w(^)/3? | hfht | 2 

SINR, =- 


(28) 
we first 


TT 


M 


1 + 


PfBi 


Hk = l,k^i Bk 

K 


Ef =1 ,^|hfh-F/3 fe 

M E£=i,jt-ii Bk 


Considering the denominator of (ITSt . using known properties 
of the inverse Wishart matrix, we have 


(29) 

We note the following properties of the terms in ( l29t . In the 
h T h* 

numerator, —> 1 as M —> oo. Also, we have 


(19) 


( 20 ) 


Hence, 

E{tr((D|H T H*D^)- 1 

where (i) = Yhf=i Pj 1 ■ Assuming that (^) converges to 
the limit, (^)°° as M,K —> oo then it can be shown that 

tr((Dg' /2 H T H*Dy 2 ) -1 ) —> (^)°°. This convergence to the 
limiting mean follows as the variance vanishes Ifl2l . Hence, 
we observe that 

«(«-!) (21) 


M 


tr(Dm r H’D^) 


—> lim 
K—*cx 

— lim 
K—> oo 


MK 
tr(D B) 
» K 
Efc=l 


j k= 1 

K 


(30) 


assuming that the limit Pk 


if the limit, f3°°, exists then in 
will also converge so that 


P°° exists. Note that 

Pk 


the terms 


K 


1 

K 


K 

E 

k=l,k^i 


Pk^P° 


(31) 


1 

Wk 


( 22 ) 


(23) 


(24) 


(25) 


Finally, from Sec. IIII-Al and (l30l >. we have 

K 

E \hfK\ 2 Pk^W- (32) 

k=l,k^i 

Combining ( l30l > and ( l32l > with ( l29l > we have the limiting SINR 
of ( |28| > for M, K —>■ oo, with fixed ^ = a, for user i given 
by 

Pfot/3? x 1 


SINR, 


=— = _ (33) 


i + smpoo xl + 

IV. Results 


where g, is the /th column of G. Under the limit operation, 
as in 0, we obtain the limit of (l25l > as M, K —»■ oo, with 


(26) 


We consider the two convergence scenarios given in Sec. 
1II-C1 and the convergence metrics given in Sec. 1III-C1 

A. Convergence Properties 

Fig. [I] shows the A ratio versus M for convergence scenario 
1 with K = 10 and K = 50. It can be seen that for both 
values of K , the A ratio will only converge to 1 when M is 
excess of 10 4 . Even for values of M as large as 100, the A 
ratio for K = 50 is more than 8 time larger than the A ratio 
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Fig. 2. MAD(E) vs K for an iid channel with a = M/K fixed 


for K = 10. However, when we consider Fig. [2] showing the 
mean absolute deviation of H T H*/M from the identity I, we 
note that this difference quickly approaches zero. For example, 
in Fig. [Q for M = 500 the A ratio is about 4 for K = 50, yet 
in Fig. [3 when K = 50, and the corresponding M = 500, the 
mean absolute deviation is less than 0.05. 

In Fig. 3, we observe that W becomes increasingly diagonal 
dominant for fixed K as M —> oo. This is because W 
has fixed dimension (K x K) and the sum of the diagonal 
elements grow faster thank the fixed number of off-diagonals. 
In contrast as both M and K grow large, W becomes less 
diagonally dominant. This follows as the number of off- 
diagonal elements increases as K 2 and the total contribution of 
the off-diagonals becomes dominant. It can be shown that the 
diagonal dominance measure grows proportionally to M 1 / 2 
for fixed K and decays proportionally to M -1 ' 2 as both M 
and K increase. For reasons of space, details are omitted. 

Considering Figs. |T| and [2] we may conclude that the mas¬ 
sive MIMO behaviour that results in a deterministic H T H*/M 
only begins to show at some very large values of M - the base 
station antenna numbers. 



Fig. 3. iid Channel Diagonal Dominance, both K fixed and = a 1 for 
M increasing 


B. Convergence Properties of ZF and MF precoders 



Fig. 4. ZF SNR and MF SINR. Single Antenna MU-MIMO, Equal Powers, 
iL = Q -l = ± 

M 10 

Figs. [3 and 0 show the convergence properties of the ZF 
and MF precoders for equal and unequal powers, respectively, 
where we plot the expected value of (IT5l) and ( 1 1 S| i for ZF 
and d25| and EH for MF, along with their corresponding 
limits. The ZF precoder in both cases of equal and unequal 
powers converges to the limit quickly. For example, in Figs. 
B] and 0 when K = 10 and M = 100 the expected value 
of the per-user SNR already approaches the asymptotic limit 
for infinite antennas. The convergence of the MF precoding 
is shown also in Figs. ED and [3 It can be observed in both 
figures that the SINR for the MF precoder is not equal to 
the SNR of the ZF precoder in the asymptotic limits; the MF 
precoder SINR is effectively reduced by a factor of Pf + 1 
when a is large (looking at (IT6t and (l26l0 . as compared to 
the ZF SNR. This difference is obvious in Figs. 13 and [3 even 
for a = 10. The MF SINR takes a longer time to converge 
because of the additional random variables in the numerator 
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Fig. 5. ZF SNR and MF SINR, Single Antenna MU-MIMO. Unequal Powers, 

— — = — 

M 10 

and denominator of (l26l >. The unequal power case for the MF 
precoder has additional terms in (l33l > which manifest itself in a 
small increase of the per-user SINR as compared to the equal 
powers. The MF precoder SINR is also less relative to ZF SNR 
due to the inter-user interference terms in the denominator of 
d23; the boost in the SINR in the numerator of (l25l) due to 
the co-phasing terms |gfg *| 2 is not enough to compensate for 
the inter-user interference given by Y2k=i,k=£i sf StSkS*- 

C. Impact of Correlation 



Fig. 6. ZF SNR and MF SINR, Single Antenna MU-MIMO, Correlated 
Channel, Equal Powers, = a -1 = 

Correlation is introduced in (0) by using the Rronecker 
model In the simulations two cases of correlation are 
shown, a high inter-element correlation of p = 0.9 and a 
low inter-element correlation of p = 0.5. Figs. [6] and [7] 
show the effect of this correlation for ZF and MF precoders 
for the mean SNR and SINR respectively. We can see that 
correlation includes a large penalty in the mean per-user SNR 
and SINR for both ZF and MF precoders (equal and unequal 
powers), when compared to the corresponding uncorrelated 



Fig. 7. ZF SNR and MF SINR, Single Antenna MU-MIMO, Correlated 
Channel, Unequal Powers, = a -1 = -k 

case. The SNR/SINRs for the lower value of correlation should 
be similar to the i.i.d. case for the two precoders - and this is 
indeed the case. 

V. Conclusion 

In this paper, we have analyzed the number of antennas 
required for massive MIMO properties. We have presented 
a method to derive ZF SNR and MF SINR for a scenario 
of unequal powers amongst the users. We have also derived 
the limit of these SNRs/SINRs as the number of base station 
antennas is increased indefinitely. We found that massive 
MIMO property that relies on small values of the off-diagonal 
elements of H T H*/M is desired, then a relatively small 
number of base station antennas will achieve convergence. 
On the other hand if the desirable matrix is to reduce the 
spread between maximum and minimum eigenvalues, then 
much larger base station antennas are needed. 

Interestingly, the per-user SNR/SINR for both ZF and MF 
precoers are less sensitive to M. In particular the per-user SNR 
for a ZF precoder converges quickly even for small values of 
M even though the SNR expression requires the computation 
of the inverse of a matrix. This leads us to conclude that for ZF 
and MF precoders, we can use the limiting value of SNR/SINR 
even for small values of M. Also, a correlated channel reduces 
MF SINR by a significant amount as compared with ZF SNR. 
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